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We present a coherent theoretical framework for computing gravitational lensing effects and redshift-space 
distortions in an inhomogeneous universe and investigate their impacts on galaxy two-point statistics. Adopt- 
ing the linearized Friedmann-Lemaitre-Robertson- Walker metric, we derive the gravitational lensing and the 
generalized Sachs-Wolfe effects that include the weak lensing distortion, magnification, and time delay effects, 
and the redshift-space distortion, Sachs-Wolfe, and integrated Sachs-Wolfe effects, respectively. Based on this 
framework, we first compute their effects on observed source fluctuations, separating them as two physically 
distinct origins: the volume effect that involves the change of volume and is always present in galaxy two-point 
statistics, and the source effect that depends on the intrinsic properties of source populations. Then we identify 
several terms that are ignored in the standard method, and we compute the observed galaxy two-point statistics, 
an ensemble average of all the combinations of the intrinsic source fluctuations and the additional contributions 
from the gravitational lensing and the generalized Sachs-Wolfe effects. This unified treatment of galaxy two- 
point statistics clarifies the relation of the gravitational lensing and the generalized Sachs-Wolfe effects to the 
metric perturbations and the underlying matter fluctuations. For near future dark energy surveys, we compute 
additional contributions to the observed galaxy two-point statistics and analyze their impact on the anisotropic 
structure. Thorough theoretical modeling of galaxy two-point statistics would be not only necessary to analyze 
precision measurements from upcoming dark energy surveys, but also provide further discriminatory power in 
understanding the underlying physical mechanisms. 

PACS numbers: 98.80.-k,98.65.-r,98.80.Jk,98.62.Py 



I. INTRODUCTION 

The standard inflationary models with a single inflaton po- 
tential predict a nearly perfect Gaussian spectrum of primor- 
dial fluctuations U |2> S B Ht] • Two-point statistics, corre- 
lation function in real space and power spectrum in Fourier 
space, constitutes a complete description of Gaussian random 
fields, and it has been widely used to understand the physics 
of the early universe from measurements of the cosmic mi- 
crowave background and large-scale structure. The recent dis- 
covery S 01 of the late time acceleration of the universe has 
spurred extensive investigations of a mysterious energy com- 
ponent with negative pressure, dubbed dark energy. Observa- 
tionally, upcoming dark energy surveys will measure galaxy 
two-point statistics with unprecedented precision from mil- 
lions of galaxies, constraining the expansion history and the 
spatial curvature of the universe. Consequently, accurate the- 
oretical modeling of galaxy two-point statistics would be cru- 
cial to take full advantage of the promise that these future sur- 
veys will deliver. 

In achieving this goal, complications arise notably from 
the nonlinear evolution of matter and scale-dependence of 
galaxy bias. In this paper we limit ourselves to the linear 
bias model [8] and study the linear theory predictions and its 
corrections, considering that recent attention has been paid 
to measu ring galaxy two-point statistics in the linear regime 
(e.g., iH 1 1 0L [ill O)- However, measurement precision is 
often highest on nonlinear scales, and proper modeling of 
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galaxy bias on nonlinear scales can substantially increase 
the leverage to constrain the underlying physics (see, e.g., 

[[IlBtHtMtBGlGl])- 

Further complication arises from the distortion of redshift- 
space structure by peculiar velocities, which results in 
anisotropy from otherwise isotropic two-point statistics Ir2(il 
Eltl . The standard practice is to analyze the angle-averaged 
correlation function or power spectrum, or to construct a lin- 
ear combination of their multipole components, suppressing 
the angular dependence of two-point statistics. However, an- 
alyzing the full anisotropic structure, though observationally 
challenging, can utilize additional information that is lost to 
some degree in the standard practice 1I22I l23l l24l l25ll . 

Gravitational lensing, often assumed to be negligible in 
galaxy two-point statistics, deflects the propagation of light 
rays, displacing the position of observed galaxies, and it al- 
ters the unit area on the sky and magnifies the observed flux, 
changing the observed number density of galaxies. The for- 
mer effect on two-point statistics is to convolve it with the 
power spectrum of the lensing potential, smoothing out the 
features in galaxy two-point statistics 02611 . The latter effect, 
known as the magnification bias ll27ll . is often used to measure 
the galaxy-matter cross-correlation function from two source 
populations separated by large line-of-sight distance 11281 129I1 . 
Recent work Il30l l3ll 13211 showed that these effects on galaxy 
two-point statistics are non-negligible at the level of accuracy 
adequate for upcoming dark energy surveys. 

However, it is unclear whether this list of additional con- 
tributions on galaxy two-point statistics is exhaustive, and 
what are the contribution terms that are ignored in the stan- 
dard method but need to be considered if higher accuracy is 
dictated by observations. Here we present a coherent theo- 
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retical framework for computing gravitational lensing effects 
and redshift-space distortions, and investigate their impacts 
on galaxy two-point statistics in an inhomogeneous universe. 
Our treatment generalizes the early work d30tl and comple- 
ments the recent work l3~ll l32tl . providing a unified descrip- 
tion of galaxy two-point statistics. However, we emphasize 
that these effects naturally arise from metric perturbations in 
our approach, comprising a complete and exhaustive set of ad- 
ditional (linear order) contributions to galaxy two-point statis- 
tic. 

The rest of this paper is organized as follows. In Sec.[TTJ we 
describe our notation for the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) metric and derive the gravitational lensing 
and the generalized Sachs-Wolfe effects. In Sec. [HI] we study 
their impacts on source galaxy fluctuations and discuss their 
correspondence to the standard redshift-space distortion and 
gravitational lensing effect. In Sec. IIV Al we derive the ob- 
served galaxy two-point statistics in real space and in Fourier 
space, and we compare the effects of each contribution term 
on the observed galaxy two-point statistics in Sec. IIV Bl We 
conclude in Sec. [V] with a discussion of the further improve- 
ment of our approach. 



n. FORMALISM 

Here we describe our notation for a background metric in 
an inhomogeneous universe and derive governing equations 
for non-relativistic matter in Sec. Ill Al Combining these with 
photon geodesic equations, we derive the generalized Sachs- 
Wolfe (Sec. Ill Bb and the gravitational lensing (Sec. Ill Q ef- 
fects, developing a coherent framework for describing how 
matter fluctuations affect observable quantities. 



for an open universe K < 0, where x is the comoving line- 
of-sight distance. The flat limit can be obtained as K — > 0. 
We will denote the covariant derivative of a three-tensor with 



,(3) 



as a vertical bar and the covariant derivative 



respect to g^ 

in the spacetime metric as a semicolon in the following. Here 
Latin indices represent 4D space-time components, and Greek 
indices run from 1 to 3, representing the spatial part of the 
metric. Throughout the paper, we set the speed of light c = 1 
We express the perturbations in the conformal Newtonian 
gauge, where ip and <fi correspond to the intuitive physical 
quantities, i.e., Newtonian potential and Newtonian curvature. 
This choice of gauge condition leaves no residual degree of 
freedom up to the first-order in perturbations. Here we only 
consider scalar perturbations, as primordial vector perturba- 
tions decay quickly in a universe with ordinary components 
and the current upper limit on tensor perturbations is order of 
magnitude smaller than the amplitude of scalar perturbations 

(e.g., JEHU!) 

Given the stress energy tensor T ab of cosmological com- 
ponents, the evolution of the matter and metric perturbations 
is governed by the Einstein equations G a b — 87rGT a b, and 
the Bianchi identities T ab = guarantee the conservation of 
energy and momentum (e.g., fl l3q,[3~7l[38l[39tl ). Current cos- 
mological observations favor a universe dominated by dark 
energy, but with non-relativistic matter as the major source 
of metric perturbations. In this universe, the scalar Einstein 
equations are 



{k 2 - 3K) 



3H 2 



a, 



3H 



(3) 
(4) 



A. Metric and Perturbations 

We assume that the homogeneous and isotropic background 
of the universe is described by the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric and its inhomogeneous 
part is represented by the perturbations for the cosmological 
fluids and the spacetime geometry: 



where 5 is the density perturbation in non-relativistic 
matter. The Hubble parameter is H — a/ a, where the overdot 
is the derivative with respect to time, dt — a dr\. The matter 
density and the Hubble parameters at the present day do are 
denoted as Q, m and Ho, respectively. The Newtonian curva- 
ture is identical to the Newtonian potential with the opposite 
sign (ijj = — (f>) in the matter-dominated era, where there 
is vanishing anisotropic stress. The conservation of energy 
momentum provides the continuity and Euler equations, 



ds 2 = -a 2 (rj) [1 + 2tp] drf + a 2 (r/) [1 + 2(f>] g ( f ) fj dx a dx (3 , 

(1) 

with the metric tensor for a three-space of constant spatial cur- 
vature K = -H 2 {1 - n ), 



r k 

d + - v 
a 



v + Hv — — ip, 
a 



(5) 
(6) 



g^dx a dx^ = d X 2 + r 2 ( X )dn 2 , (2) 

where a(rj) is the scale factor for the expansion of the back- 
ground as a function of the conformal time r\, and the comov- 
ing angular diameter distance is r(x) — K~ x / 2 s\y\{\/Kx) 
for a closed universe K > and (— A') -1 ' 2 sinh(V— K\) 



where v is the velocity of non-relativistic matter in units of c. 
In the conformal Newtonian gauge, the relativistic equations 
on sub-horizon scales correspond to the usual Newtonian 
equations, 

VV = ^ tl m 6 -, (7) 
2 a 
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v = — 



2 

3 fl 



a 2 Hf 



Vfl>, 



(8) 



and perturbations in the geodesic equations for the temporal 
and spatial components are 



where / = din D /din a and D is a growth factor of 
the matter density perturbation. The evolution of the density 
perturbation is related to the Newtonian potential by 



8 + 2H8 = V>, 



(9) 



and the growth factor D is a growing solution of this differ- 
ential equation, normalized to a unity at oq. Full relativistic 
consideration results in additional multiple terms in the right- 
hand side of the equation (e.g., iH), but they are suppressed 
at least by the ratio of a characteristic scale 1/fc to the Hubble 
distance 1/H. Note that we have interchangeably expressed 
equations in Fourier space and configuration space, which is 
valid to the linear order in perturbations and significantly sim- 
plifies the manipulations. 



B. Geodesic Equations and Sachs- Wolfe Effects 

The propagation of light rays is described by a photon 
geodesic x a (X) with an affine parameter A, and a null vector 
k a = dx a /d\ tangent to x a is determined by the null equation 
(d.s 2 = k a k a = 0) and the geodesic equations (k a . k b = 0). 
In a perturbed FLRW universe, the null vector can be ex- 
pressed as 



k" = - (1 +8v), k a 
a 



{e a + Se a ), (10) 



where v and e a are the photon frequency and its (time- 
reversed) propagation direction from the observer, and the di- 
mensionless quantities 8v and 8e a represent their perturba- 
tions. In a homogeneous expanding universe, the null vec- 
tor follows the usual relations v oc 1/a, e a e a = 1, and 
de a jdr\ = e$e a \p, and indeed Eq. ( fTOb may be derived from 
the null and the geodesic equations. For a comoving observer 
whose rest frame has vanishing energy flux, the four velocity 
is u a = (1/a, 0) and the observed frequency v a \ }S of a pho- 
ton source is related to the frequency v c at the emission by a 
redshift parameter, 



1 + z 



(k a Ug) e _ _Uc_ 
(k a U a ) 

obs ^obs 



(ii) 



where we assumed a Q b s = an = 1. 

In an inhomogeneous universe, the observed redshift £ b s 
deviates from the true redshift z. Perturbations in the null 
equation is 



dy drj 



(13) 



dy 



(8e a +2<t>e a ) = 8e p e c 



de a 



where we used the zeroth order null geodesic 
d/dy = d v — e a d a = (a/u)(d/dX) and kept the terms 
to the first order in perturbations. 

The four velocity of a comoving observer is now u a — 
((1 — ip)/a, v a /a) and the observed redshift is 



(k a u a ) c 



(1 + z) [l + {8v + i, + v a e a f a ]. 

(15) 



(k a U a )obs 

This can be further simplified by using Eq. ( fT3l ) as 



1 + z ohs = (1 + z) x 



l + V(z)~V{0) 



(16) 



iP(z) + yj(0)+ I" dy-^-(cP-iP) 
Jo 9-q 



where V — v a e a is the line-of-sight velocity lf30l l40l |4lll . 
The additional terms in the square bracket alter the simple 
redshift-distance relation in Eq. dTTb . giving rise to the 
standard redshift-space distortion by peculiar velocities, 
the Sachs-Wolfe effect by gravitational redshift, and the 
integrated Sachs-Wolfe effect by the time evolution of gravi- 
tational potential across which photons propagate. Hereafter 
we will collectively refer to these effects as the generalized 
Sachs-Wolfe effect. 



C. Gravitational Lensing 

In a homogeneous universe, the gravitational lensing effects 
vanish and light rays propagate with the direction unchanged. 
For a photon source at z-axis in an inhomogeneous universe, 
the propagation direction from the observer is n = e a = 
(0, 0, 1) and the null vector is k x ^ = -(u/a)8e x ^. The null 
vector is further related to the photon position r(x)h — (x, y) 
on the sky at any time by 



k x, y _ A( rfi ) = -^-(rn), 
d\ a dy 



(17) 



e a 8e a = 8v + ip — <f>, 



(12) 



where we replace the derivative with respect to the affine pa- 
rameter by using the zeroth order null geodesic, consistent to 
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the first order in perturbations. The spatial component of the 
geodesic equation (Eq. |[T4"l ) is then 



^{8e x -y) = --^r(rn) - V> - 4> W = 2 ^ a ■ 
ay ay z 



(18) 



Since gravitational lensing conserves the surface bright- 
ness, the observed surface brightness J b s (n) on the sky is 
simply the intrinsic surface brightness at the source position 
s: iobs(n) = I(s), and the source position § can be obtained 
by integrating Eq. ( fT8l along the photon geodesic 



s = n + V *(n), 



with the projected lensing potential 



m = -2[ y "dy' [ V dy-^L-^ 
Jo Jo r (Xs) r(x) 



(19) 



2 r dy ^(y) 7 ) X \ ?\ , (20) 
Jo 



r(Xs) r(x) ' 



where V is the derivative with respect to n, and 
X s = Jq 3 dz/H(z) is the comoving line-of-sight dis- 
tance to the source redshift z s . The integration along the 
unperturbed photon geodesic dy is often called the Born ap- 
proximation. Following the literature, we take the geodesic as 
the photon radial direction d\, but note that d/d\ = d v — d x . 

The convergence ft(n) is defined as V 2 ^>(n) — — 2«(n) 
and it is further related to density fluctuations along the 
geodesic by Poisson's equation (Eq.JT)) 



/c(n) 



d X (V 2 



V 2 ) ip[r(x)n,x] 



r(Xs ~ x) r(x) 



Xs dx S[r(x)n,x] r(x 

a(x) r( X 



r (Xs) 

X) r(x) 



(21) 



The contribution from the radial derivatives V 2 , is pro- 
portional to the potential difference between the source and 
observer, and this boundary term is negligible compared to 
the first term ll42l l43ll . Numerical ray tracing experiments 
through TV-body simulations show that the weak lensing 
approximation to the first order in perturbations is accurate 
even in nonlinear regime when nonlinear matter power 
spectrum is used in place of linear matter power spectrum 
II42I1 . Also note that all the prior results for a single source 
redshift can be readily generalized to a source population 
with a redshift distribution W(xs) by integrating the results 
over Xs with W(xs) in the integrand. 

While conservation of surface brightness guarantees that 
photons are neither destroyed nor created, gravitational de- 
flection distorts the cross-section of a bundle of light rays, 



magnifying (or de-magnifying) observed fluxes. Gravitational 
lensing magnification /i(n) is related to the Jacobian of a map- 
ping from the image plane to the source plane by 



d 2 s 



d 2 h 



I + VVf(n; 



(22) 



[l- K (fi)] 2 - 7 2 (n) , 



where I is a unit 2x2 matrix and 7(n) is the tangen- 
tial shear. In the weak lensing regime, /i(n) = 1 + 2 «(n). 

Gravitational lensing also modifies the propagation time of 
light rays in two ways, compared to the light travel time in the 
absence of the gravitational lensing effects: it distorts the pho- 
ton geodesic, increasing the path length that photons travel, 
and the gravitational potential retards the light travel time. 
The former is referred to as the geometric time delay [44] 



/.x 1 r(xi) r(x 
T-geo(n) = - 



2 r{ Xs 



XI) 



V*(n)-W(n), (23) 



and the latter is the potential or Shapiro time delay [45 



r{Xs - Xi) 



(24) 



These effects can be derived by using the small angle approxi- 
mation in deflection and the relation dr\ = (1 — tp+ <fi)dx from 
the metric in Eq. ([T}. Note that the proper time delay can be 
obtained by multiplying the lens redshift 1 + z\ in the limit of a 
single lens case, and this derivation in a cosmological context 
recovers the standard relation for time delay. 



III. SOURCE FLUCTUATIONS 

Inhomogeneous matter fluctuations in the universe deflect 
the propagation of light rays, giving rise to the gravitational 
lensing effects. The generalized Sachs-Wolfe effect also arises 
from the same matter fluctuations responsible for the gravita- 
tional lensing effects. Having discussed the basic mechanism 
of the gravitational lensing and the generalized Sachs-Wolfe 
effects that complicate the simple interpretation of observable 
quantities, we now investigate their impact on an observed 
overdensity field £ b s (n, z) of source galaxies. Contributions 
to <5 bs(n, z) come from matter fluctuations in addition to the 
intrinsic overdensity S(h, z) of source galaxies. Noting that 
the contributions can be linearized and added to the first or- 
der in perturbations, we separate these contributions as two 
physically distinct parts: one that involves the change of vol- 
ume, and one that involves the intrinsic properties of source 
galaxies. The impact on galaxy two-point statistics will be 
discussed in the following section. 
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A. Volume Effect 

Consider a unit comoving volume dV = r 2 (x)d^ldz / H (z) 
and a unit flux interval df, and let n(n, z, f) be the comoving 
number density of source galaxies. The generalized Sachs- 
Wolfe effect alters the unit comoving volume dV. Note, how- 
ever, that it not only changes the unit redshift interval dz, but 
also changes both the angular diameter distance r(x) and the 
Hubble parameter H(z). By imposing the number conserva- 
tion, the observed number density of the source galaxies can 
be obtained by 



«-obs(.2obs) = n(z) [1 + SV] 



r 2 (x) H(z ohs ) dz 
r 2 {Xobs) H(z) dz ohs : 



(25) 



where 5V — (2<j) + e)% represents the distortion of volume el- 
ement, when it is transformed from the conformal Newtonian 
gauge to the local Lorentz frame, where the velocity of non- 
relativistic matter vanishes. We give a more rigorous deriva- 
tion in Appendix [B] The solid angle dfl remains unaffected 
by the generalized Sachs-Wolfe effect. 

If the mean comoving number density evolves slowly com- 
pared to the redshift change due to the generalized Sachs- 
Wolfe effect n(z) — n(z b s ), contributions to <5 bs(2 bs) arise 
solely from the change in volume element dV, 



<5obs(zobs) = S(z)- 



-SV, 
(26) 



where we rewrote Eq. ( fTBT ) as 1 + z b s = (1 + z)(l + e) and 
the contribution e from the generalized Sachs-Wolfe effect is 



if we ignore the Sachs-Wolfe and the integrated Sachs-Wolfe 
effects in Eq. ( f27b . This recovers the standard relation for 
redshift-space distortions H S 53] • Note that the stan- 
dard method ignores the contributions in Eq. (Eot from the 
Sachs-Wolfe and the integrated Sachs-Wolfe effects. We dis- 
cuss their impact in Sec.lIVI 

Gravitational lensing magnification increases the flux inter- 
val df and the solid angle dft by a factor of \i, respectively. 
With the number conservation in dV and df, the observed 
number density is therefore 



n bs(/obs) = n(f) 



df dn 



df ohs dCl 



obs 



(29) 



Similarly, if the mean comoving number density is the same 
over the flux change due to lensing magnification (i.e., the 
source luminosity function is flat), the observed overdensity 
is then 



(5 bs(n) = 5(n) - 4 n(n), 



(30) 



reflecting the change in volume and flux. 

Gravitational lensing displaces the source position on the 
sky according to Eq. ( fT9l , and the observed number density 

is n bs(n) = n n + V<f(n) . By Taylor expanding n Q b s (n) 

to the first order in ^(n), the observed overdensity can be 
written as 



Wn) = 5(n) + W(n)-V5(n). 



(31) 



s(z) = V(z) - V(0) - ip(z) + V>(0) 



x , dtp 
or] 



(27) 



In the Einstein-de Sitter universe, the Newtonian potential 
is constant and hence the integrated Sachs-Wolfe effect van- 
ishes. In general, as we show in the next section, the peculiar 
velocity effect is dominant over the Sachs-Wolfe and the inte- 
grated Sachs-Wolfe effects, and e(z) ~ V{z) - V(0). Note 
that while our derivation so far is valid for nonflat universes, 
in deriving Eq. d26*l i we assumed that the spatial curvature K is 
close to zero. The second term in Eq. ( f26b has a multiplicative 
factor VKx/ tan(VXx) f° r a closed universe K > and 
V 1 —Kxl tanJi(v— ~Kx) f° r an open universe K < 0, which 
becomes a unity as K — > 0. 

With a proper line-of-sight distance r p = x{ z ) /(I + z ) an d 
a normalized peculiar velocity u = V(z)/ H(z), Eq. d26b can 
be rearranged as 



<Sobs(Zobs) = S(z) 



2u 



du 
dr p ' 



(28) 



Note that the additional contribution is already in the second 
order in perturbations and furthermore it vanishes on average, 
because the deflection angle V'I'(n) has no preferred direc- 
tion. The first non-vanishing effect from gravitational lensing 
displacement comes in the second order in ^(n) 13111 . and we 
therefore ignore this effect. 

Finally the gravitational time delay decreases the arrival 
time of photons in an overdense region, compared to that in 
the absence of lensing. The net effect is therefore that we 
sample sources at farther distance in the fixed time interval 
ll48ll . However, for discrete sources the effect vanishes as long 
as the life time of the sources is longer than the time delay. 



B. Source Effect 

The generalized Sachs-Wolfe and the gravitational lensing 
effects modify a unit volume and a unit flux interval, lead- 
ing to the contributions to £ b s (n, z, /). Furthermore, the 
changes in observed redshift and flux can result in differ- 
ent mean number densities, if the redshift distribution of the 
source galaxy population varies in the redshift interval or the 
luminosity function is non-trivial over the flux change. These 
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additional contributions from the change in mean number den- 
sities are related to the intrinsic properties of source galaxies, 
and we collectively refer to these effects as the source effect. 
However, note that while the source effect may be absent for 
some galaxy populations, the volume effect is always present. 
Therefore, we keep together the contributions from the vol- 
ume effect in considering the source effect. 

We first consider the effect of gravitational lensing magnifi- 
cation. Lensing magnification not only increases dil and df in 
Eq. (|29] l, but also changes the number count of source galax- 
ies, if the luminosity function is non-flat, i.e., n bs(/obs) 7^ 
n(f). Assuming n(f)df oc f~ s df with a constant slope s 
over a narrow flux range df, the observed number density can 
be expressed as 



_ v n(f obs /fi) s _ 2 
n bs(/obs) = 5 = n W°t,s) A* 



(32) 



and the observed overdensity is now 



<Sobs(n) 



S(h) 
5(h) 



(2s- 
5(p- 



4)«(n) 
0.4)«(n) 



(33) 



where we used the logarithmic slope p = d\ogn(m)/dm — 
0.4(s — 1) in a sample with limiting magnitude to. In the 
literature, these contributions from both the volume and 
the source effects are referred to as the magnification bias 

113, miUHim]. N ° te ^ ^ biaS Can be either P ositive 
or negative, depending on the slope p, and the volume effect 

can be canceled by the source effect with p — 0.4 (see lEill 
for the recent detection from the Sloan Digital Sky Survey). 

The redshift distribution of source galaxies also affects 
the mean number counts due to the generalized Sachs- 
Wolfe effect. For a redshift distribution n(z)dz oc 
z a exp [— (z/zq)P~\ dz, the observed overdensity can be ob- 
tained by substituting n(z) with n [z ohs - (1 + 

^obs ) £ J 1 



<5obs(z) = S(z) 



l + z 



(34) 



where the second term in the right-hand side is the ad- 
ditional contribution related to the evolution of source 
galaxies, and the rest of the additional terms come from the 
volume effect in Eq. 



We then compare their impact on galaxy two-point statistics 
in Sec. EV1 

For a sample of galaxies at redshift z selected with a lim- 
iting flux / and narrow intervals of dz and df, the observed 
overdensity <5 b s (n, z, f) is the sum of the intrinsic overden- 
sity field 5(h, z, f) and the contributions from the gravita- 
tional lensing and the generalized Sachs-Wolfe effects: 



5 obs (n, z, /) = 6 + 6 mh + 6 Z + S e 



From Eq. (l33l , the magnification bias is defined as 



(35) 



6 mh (n, z,f) = 5 [ P (f) - 0.4] K (fi, z), (36) 



with redshift z being the source redshift of the convergence 
«(n) in Eq. (EH . Considering e(z) ~ V(z) — V(0), we call 
the volume effect in Eq. ( |26] i as the redshift-space distortion 
bias, 



S x (n,z) = -2 



= -2 



1 



l + z , TT d/e 

-W £ - {1 + z)H dz{H 

l+z l+z dH 

e H e 

H X H dz 

- z de CTr 
£ + SV. 



-e + SV 



H d X 



(37) 



Note that the generalized Sachs-Wolfe effect e(z) im- 
plicitly depends on the direction fi via the line-of-sight 
velocity V(z) — v a e a — h ■ v(n, z), but it is independent of 
the limiting flux /, provided that galaxies have no velocity 
bias (i.e., galaxies and matter follow the same velocity field). 
Finally, the evolution bias is defined from Eq. d34b as 



6 evo (n,z,f) = - 



l + z 



(38) 



where the directional dependence comes from e and the evo- 
lution coefficients (a, f3, Zq) depend on the galaxy sample se- 
lected with the limiting flux /. While the evolution bias aris- 
ing from the difference between n(z) and n(z b s ) was rec- 
ognized i20l I22I l47ll . it has been ignored in the literature. 
However, we show in Sec. [IV] that the evolution bias can 
be significantly enhanced. Last, we want to emphasize that 
equation (l35l l is gauge-invariant as is written in the conformal 
Newtonian gauge. 



C. Summary 

We have investigated the effects of inhomogeneous matter 
fluctuations on observed overdensity fields. Here we summa- 
rize their contributions and clarify the functional dependence. 



IV. GALAXY TWO-POINT STATISTICS 

We have derived additional contributions of the gravita- 
tional lensing and the generalized Sachs-Wolfe effects to 
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the intrinsic density fluctuations in Sec. [Ill] fully consis- 
tent up to the first order in perturbations. Given two sam- 
ples of galaxies with limiting fluxes /i and f%, the observed 
galaxy correlation function is then £ bs(ni, Z\, tl 2 , z 2 ) — 
(<5 bs(ni, Z\) <5 bs(n2, Z 2 )) and the observed power spectrum 
is (<5 bs(ki,zi) <5* bs (k 2 ,z 2 )) = (2vr) 3 ^(ki - k 2 )P obs (fci). 
In Sec. lIV Al we derive this ensemble average of all the com- 
binations of each component in <5 bs in Eq. ( T35l l. after we sim- 
plify the equation. We then discuss their impact on the ob- 
served galaxy two-point statistics by analyzing specific exam- 
ples in Sec. UVBl 



than 8k by the ratio of correlation scale 1 /k to the Hubble dis- 
tance 1 /H or the line-of-sight distance x (roughly of order 
1 / H), except the third component: the partial derivative with 
respect to \ cancels the correlation scale 1/k and hence the 
amplitude of the third component is of order 6k, larger than the 
other components in the redshift-space distortion bias. There- 
fore, we only keep the dominant component in the redshift- 
space distortion bias 13011 , 

xr* ^ 1 + zdV 
5 z (n,z)c,-—-, (41) 



A. Correlation Function and Power Spectrum 

Here we compute the observed galaxy correlation function 
£obs and power spectrum P oos - However, as some compo- 
nents in S bs are smaller than other components, their combi- 
nations are even smaller by an order-of-magnitude. We there- 
fore start by estimating the auto-correlation functions of each 
component and simplify the equation before we compute all 
the cross-correlation functions and power spectra. 

We first consider the correlation of the redshift-space dis- 
tortion bias £ zz = (S z (hi,zi) <5 z (n.2, #2))- The redshift-space 
distortion bias S z in Eq. d37l i has five components that depend 
either e(z) or its partial derivative with respect to z or %, and 
the contribution e(z) from the generalized Sachs-Wolfe effect 
in Eq. d27| i has also three different components that depend 
on the peculiar velocity, the Newtonian potential, and its time 
derivative. The Newtonian potential and the peculiar velocity 
in Eqs. (Q and © take the simple form in Fourier space 



^k 

Vk 



2 

iH fa (5k 



_ 4 

a k 2 ' 
k 



(39) 
(40) 



On a typical correlation scale 1/k, they scale as H 2 Sk/k 2 
and HSk/k with / ~ 1 at z ^ 1: ipk is smaller than Vk by the 
ratio of the correlation scale 1/k to the Hubble distance 1/H. 
Similarly, the integrated Sachs-Wolfe effect is of the same 
order as the Newtonian potential and it vanishes in the limit of 
zero cosmological constant, i.e., Einstein-de Sitter universe, 
because it is proportional to the time derivative of the ratio 
of the growth factor to the expansion scale factor D(z)/a. 
Therefore, we can safely ignore the Sachs-Wolfe and the 
integrated Sachs-Wolfe effects and we assume e(z) ~ V(z). 
Note that given a particular realization of the observer's rest 
frame, its peculiar velocity V^O) is uncorrelated and the 
unobservable potential t/j(0) in Eq. ( 1271 ) can be absorbed by a 
gauge transformation. 

With the assumption e(z) ~ V(z), we further simplify 
Eq. d37l i by comparing the five components in the redshift- 
space distortion bias, and similar justification was made in 
IBOIl . Respectively, each component scales as 5k /k\, HSk/k, 
dSk/kd\, H5k/k, and H5k/k, and hence they are smaller 



consistent with the standard relation for the redshift-space dis- 
tortion, justifying its nomenclature. However, note that all 
these ignored components are proportional to e. At low red- 
shift, they contribute to galaxy two-point statistics at the sub- 
percent level, while we show in Sec. IIV Bl that at higher red- 
shift their contribution is somewhat larger. 

Having substantially reduced the number of combinations 
for an ensemble average, we are now well positioned to com- 
pute correlation functions and their power spectra. For two 
galaxy positions xi = [r(xi)ni, %i] andx 2 = [r(x 2 )n 2 , x 2 ], 
the auto-correlation of the redshift-space distortion bias is 



(4(ni,2i) S z [h 2 ,z 2 )) (42) 
d 3 k 



/l/2 /(^ e " (X " X2)Pm(fc;Zl ' Z2) F 
r°° dfc ^3 

hh L i^ a;zi,Z2) 

gjo(fcr 3 )P (7) ~ ^ 2 (fcr 3 )P 2 ( 7 ) + ^(kr^P^) 

where the 3D comoving separation is r 3 = 

[r(x) 2 A0 2 + ( X2 -Xi) 2 ] 1/2 with Ad = |nj - n 2 | and 
X — (xi + X2V2, and the angle subtended by the comoving 
separation is 7 = cos A6 = (X2 - Xi)/r 3 - p n(x) and j n (x) 
are the n-th order Legendre polynomial and spherical Bessel 
function, respectively. We assumed the distant observer 
approximation such that k z is the line-of-sight component 
of the wavenumber k, but it can be relaxed by replacing fc* 
by [(rii • k)(n 2 • k)] . The linear matter power spectrum is 
computed by P m (k; Z\, 22) = D(z\)D(z-2)P m (k), while 
we use P m (k; zi, Z2) — P m (k;z) with z = (z\ + z 2 )/2 
when we compute the effect of the nonlinear matter power 
spectrum using the Smith et al. 115 311 approximation. The 
power spectrum of the redshift-space distortion bias can be 
readily read off from Eq. d42l and its power is boosted along 
the line-of-sight by /i/ 2 /i^ with Hk — k z /k. 

Next we consider the correlation of the evolution bias. The 
observed redshift z oos is different from the true redshift z due 
to the generalized Sachs-Wolfe effect and the redshift distri- 
bution of the source mean number density gives rise to the 
evolution bias. The evolution bias S cvo is proportional to 
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e(z) ~ V(z) and it is typically smaller than S z by the ra- 
tio of a correlation scale 1/fc to the Hubble distance 1/H. 
However, beyond the mean redshift of source populations, the 
mean number density changes exponentially and the evolution 
bias can be substantially boosted by the prefactor 



where J n (x) is the n-th order Bessel function of the 
first kind. Assuming that the source redshifts are suf- 
ficiently high and hence £ m b is independent of z% and 
z 2 , the power spectrum of the magnification bias is 

P mb = (2n)S D (k z )(5 Pl - 2)(5p 2 - 2)r 2 (x)C^ k±r{& 
where the angular power spectrum of the convergence is 



E(z;f) = - 



l + z 



a-(3[- 



(43) 



fir. 



dx 



r(x - x) 



1 2 



o(x) r (x) 



Pn 



k = 



I 



defined such that Eq. d3~8l becomes <5 0VO = E(z) e(z). While 
the exact functional form of E(z) depends on the assumed 
redshift distribution, it captures the general trend of the en- 
hancement in S cvo beyond the mean redshift. The correlation 
of the evolution bias is therefore 



r(x) ' 
(47) 



£cvo — \0cvo(fil> z l) ^evo 

'd 3 k 



(44) 

2 



(HfaEUHfaEh LiJl e *-C«x-*») P m (k; Zi,z 2 )^ 
f°° dk k 

(HfaEUHfaEh p m (k; Zl ,z 2 ) 

Jq k Ztt 



^o(fcr 3 )Po(7)-^2(fcr 3 )P 2 (7) 



where the subscripts in the round brackets represent that 
the products (HfaE) are computed at z\ and z 2 . Its power 
spectrum is also anisotropic and has structure similar to the 
redshift-space distortion bias. 

Finally, the inhomogeneous matter fluctuations along the 
two lines-of-sight result in the correlation of the magnification 
bias 



The Dirac delta function results from our assumption that 
£ m b is a function of transverse direction only, but it can be 
somewhat relaxed by replacing (2ir)5 D (k z ) by a survey win- 
dow function 05511 . Note that while we are interested in how 
the magnification bias affects the 3D correlation of the in- 
trinsic source fluctuations, the magnification bias arises from 
the matter fluctuations along the line-of-sight (not at a sin- 
gle redshift plane) and thereby angular correlation function 
and its angular power spectrum are better suited for quan- 
tifying its statistics. Indeed, the correlation function of the 
magnification bias is identical to the angular correlation func- 
tion, £ m b(ni, zi, n 2 , z 2 ) = w mh (A6; z x , z 2 ), and we relate 
2D angular power spectrum to 3D power spectrum by P(k) = 
(2ir)6 D (k z )r 2 (x)Ci =k±r{x) (see Appendix©. 

With all the additional contributions of the gravitational 
lensing and the generalized Sachs-Wolfe effects in hand, the 
correlation of the intrinsic fluctuation of sources is modeled 
using the linear bias model, 



&nt = (<5(ni,zi) 8(& 2) z 2 )) 



(48) 



h b 2 



£mb = (<Jmb(ni,2l) S mh (h 2 ,Z 2 )) 

= (5pi-2)(5p2-2 



(2tt)3 
dk k 3 



D ik-(xi-x 2 ) 



Pm {k 

; zi,z 2 ) 



V 2 

r[xi ~ X) r(x2 ~ x) 



dx 



r(x) 



(45) 

2 



a(x) 



where we used the Limber approximation 115411 (see Ap- 
pendix©. Without loss of generality, we assumed Z\ < z 2 . 
The projected correlation function w p (R) is obtained by 
integrating the 3D matter correlation function £ m (x) along 
the line-of-sight at a fixed redshift z(x) and 2D transverse 
separation R, 



f dk k 

hhj —— P m (k;z l7 z 2 ) jo(kr 3 ), 



where the constant linear bias factors b\ and b 2 are the 
ratio of the intrinsic source fluctuation to the underlying 
matter fluctuation at Z\ and z 2 . 

To complete our calculations of £ b s , we now compute the 
cross-correlation functions and power spectra between the in- 
trinsic source fluctuation and the fluctuations from the gravita- 
tional lensing and the generalized Sachs-Wolfe effects. First, 
the redshift-space distortion bias and the intrinsic source fluc- 
tuation provide two cross-correlation functions 



u p [R; z] = 



dr\\ £ 

k dk 
~2^~ 



r 3 = ijR 2 + rf,z 
P m (k;z)J {kR), 



(46) 



£ fe = (<S(rii,zi) <5 z (n 2 ,z 2 )) 
d 3 k 



(49) 



= hh 
= hf 2 



k-(xi-x 2 ) p (u. \ _z_ 

(27r)3 . Pm{k,Zx,Z 2 ) fe2 

50 dk k 3 



k 2tt 2 



P m {k;z\,z 2 ) 



~jo(kr 3 )P (~/)-lj 2 (kr 3 )P 2 ( 7 ) 
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and similarly for £ z s — (<5 z (ni,zi) #(112,22)) with the 
two indices exchanged in Eq. d49i >. Combined with £ zz in 
Eq. d42l >. these two cross-correlation functions constitute the 
standard redshift-space correlation function 

£z-dist = £int +6» + = (50) 

2=0.2,4 



which is often expressed in terms of the multipole components 



r°° dk k 3 

& = ciiPufh) h b 2 i / — —2 P m (k;zi,z 2 ) ji(kr 3 ), 

(51) 



redshift, i.e., £ m b ,5 = Q. The power spectrum is also related 
to the angular power spectrum of the cross-term 



C, 



8 mb 



61 (5 P2 - 2) [ 



(56) 



v r(X2 - Xi) p 
on n r 2 



as F,5 mb = (2n)6 D (k z )rjCf = f^ ri . Finally, all the 
cross-correlations that involve S cvo are zero, since S cvo is odd 
in the line-of-sight component V of peculiar velocities and 
the universe has no preferred direction. Two remaining cross 
terms £ z m b and £ m b z also vanish, since S z is proportional to 
k 2 z and the line-of-sight fluctuations are smoothed out in the 
Limber approximation. 



with its coefficients 



B. Comparison 









■I 


w 





1 , 01+02 

1 T 3 



Pi 02 



\{01 + &) + f 01 02 



(52) 



8 0iP2 



where (3 = f /b. Analogously, the redshift-space power spec- 
trum is 



P 



z— dist 



= Pint + PzS + Psz + Pzz (53) 
= [l + (01+0 2 )»l+010 2 4]Pint(k) 

= Pi(»k)Pi(k), 



(=0,2,4 



with the intrinsic source power spectrum P; n t(fc) = 
bi b 2 Pm(k; Zi, z 2 ), and its multipole components Pf(k) in 
Fourier space are related to the multipole components £i(rs) 
in real space as 



/>oo 

Pf(k) = 4%i l / dx x 2 ii{x) ji(kx). (54) 
Jo 



Since the magnification bias arises from the matter fluctu- 
ations along the line-of-sight, it correlates with the intrinsic 
source fluctuation at lower redshift (z\ < z 2 ), 



Umb = (S(ni,zi) 5 mh (h 2 ,z 2 )} 

= bi(bp 2 - 2) —^-n. m w p 7-1 A0; 

V 2 I air 2 



but the correlation vanishes when the source is at higher 



To compare the additional contributions to the observed 
correlation function £ b s and power spectrum P b s , we con- 
sider near-future dark energy surveys that will target galaxies 
and quasars to measure their correlation function and power 
spectrum at high redshifts. For example, the baryonic os- 
cillation spectroscopic survey (BOSS) will measure 1.5 mil- 
lion luminous red galaxies to determine the angular diame- 
ter distances at z = 0.35 and 0.6, and use 160,000 quasars 
to measure the clustering of Lyman-a forests at z = 2.5 
H571 15811 . For the purpose of illustration, we show our calcu- 
lations of galaxy two-point statistics at these redshifts. Here 
we adopt a flat ACDM universe: the cosmological parame- 
ters are the matter density Vl m h 2 = 0.137, the baryon density 
n b h 2 = 0.0227, the Hubble constant h = 0.70, the spectral 
index n s = 0.96, the optical depth to the last scattering sur- 
face t = 0.084, and the primordial curvature perturbation am- 
plitude A 2 = 2.457 x 10" 9 at k = 0.002 Mpc" 1 (correspond- 
ing to the matter power spectrum normalization a% — 0.817), 
consistent with the recent results (e.g., lf33i[34l[35Tl ). The mat- 
ter transfer function is computed by using cmbfast |[59tl . 

Figure Q] examines the separate contributions of the grav- 
itational lensing and the generalized Sachs-Wolfe effects to 
the observed two-point correlation function of galaxies. We 
show the correlation functions of the intrinsic galaxy fluctu- 
ations (£i n t/6 2 ; solid) and the redshift-space distortion bias 
(£zz; dotted), and their cross-correlation function (£,sz/b = 
£,zs/b; long dashed). The correlation functions are computed 
by using the linear (thin) and the nonlinear (thick) matter 
power spectrum. Note that they only differ on small scales 
and the nonlinear effect decreases at high redshift as shown in 
(5<j-)Fig.[T]2 to Fig.QJi, going from z = 0.35 to z = 2.5. The source 
galaxies are assumed to be at the same redshift (z = z\ = z 2 ) 
Z A shown in the figure legend, and thus 3D separation r 3 is equal 
to 2D projected separation R = r(x)A0. However, two 
galaxy populations are separately placed at zi = 0.35 and 
z 2 = 0.6 in Fig.[T]i, and the x-axis represents projected sepa- 
ration R, rather than 3D separation 7-3. 
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separation (h 'Mpc) projected separation (h 'Mpc) 



FIG. 1: Dissection of the observed two-point correlation function of galaxies. Solid, dotted, and long dashed lines represent correlation 
functions of the intrinsic source fluctuations £i n t/& 2 and the redshift-space distortion bias £ zz , and their cross-correlation function £gz/b = 
£z.s/b, respectively. Correlation functions of the magnification bias £ m b/(5p — 2) 2 and the evolution bias ^ cvo /E 2 are shown as short dashed 
and short dot-dashed lines. Note that while the galaxy bias factor b and magnification bias factor (5p — 2) are of order unity, the evolution boost 
factor E can be an order of magnitude larger (see Fig. [3}- The correlation functions are computed by using the linear {thin) and the nonlinear 
{thick) matter power spectra, and source galaxies are assumed to be at the same redshift indicated in the legend. The correlation functions of 
the intrinsic source fluctuations become negative at r% > 128/i _1 Mpc, where its absolute value is plotted. Panel {d) plots the cross-correlation 
function £g m b/foi(5p2 — 2) of the intrinsic fluctuation of source galaxies at Z\ = 0.35 and the magnification bias from source galaxies at 
Z2 = 0.6 as long dot-dashed lines. With large line-of-sight separation 600/i _1 Mpc, only £ m b and £5 m b that depend on projected separation 
R rather than 3D separation r$ itself are appreciable, i.e., £j n t ~ (,sz — £z6 — £zz — £cvo ~ 0. 



The solid lines £;„t / b 2 are identical to the matter correlation 
function £ m and the linear bias factor b is constant. How- 
ever, the nonlinear evolution and galaxy formation process 
complicate the relation between galaxies and underlying mat- 
ter fluctuations, and galaxy bias becomes scale-dependent on 
small scales, even when the nonlinear matter power spectrum 
is used (e.g., Il60l I61I0 . While we plot the correlation func- 
tions at r3 ~ 0.5 — 200/i~ 1 Mpc for completeness, the valid- 
ity of our calculation is limited to the linear regime. The solid 
lines at r$ = 151 Mpc (= 106/i _1 Mpc) show prominent en- 
hancement in the clustering amplitude, known as the baryonic 
acoustic peak ll62j|63ll . The baryon-photon plasma in the early 
universe propagates as sound waves and these periodic oscil- 
lations in Fourier space translate into one peak in real space 
with its width deviating from a sharp delta function due to the 



termination of the harmonic series, determined by the horizon 
size at the cosmological recombination epoch. Note that the 
correlation function becomes negative at r$ ~ 128/i _1 Mpc, 
beyond which we plot its absolute value. 

The correlation functions (£ zz , £$ z , and £ z g) of the redshift- 
space distortion bias have the overall shape similar to £i„t- 
However, since £ zz , ^ z , and £ z ,5 in Eqs. d42b and d49l ) have ad- 
ditional functional dependence on spherical Bessel functions 
j2(x) and ji(x) compared to £i nt in Eq. 148K it puts more 
weight on higher k and hence the nonlinear effects persist 
up to r3 ~ 10/i -1 Mpc in Fig.Q]}, larger than 3/i _1 Mpc for 
£i n t. However, the incoherent superposition of the additional 
Bessel functions washes out the acoustic peak in the correla- 
tion functions of the redshift-space distortion bias, leaving lit- 
tle structure in £ zz , £a z , and £ z s at the acoustic scale. Since the 
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FIG. 2: Multipole components of the redshift-space correlation func- 
tion £ z _dist and the evolution bias £ cvo at z = 0.35. We define multi- 
pole components of £ e vo as £ cvo = £° vo (?"3 ) Po (7) + £cvo (^3) P2 (7), 
in the same way multipole components of £ z _dist is defined in 
Eq. d50t . We assume that the galaxy bias factor is b = 2 and the evo- 
lution boost factor is E = 100 for a proper comparison. The corre- 
lation functions are computed by using the linear matter power spec- 
trum only. £o(r-i) becomes negative at rs > 128fo _1 Mpc, where its 
absolute value is plotted. The inset shows the correlation £i nt of the 
intrinsic source fluctuations around the acoustic scale and it is related 
to the monopole by £0 = £lnt ■ c o in Eq. l |51t . where Co = 1.24 in 
our fiducial model. 

observed correlation function £ Q b s is the sum of all the contri- 
butions and it is hard in practice to separate each contribution 
from one another, it may look as if £; nt is swamped by £ zz , 
£i z and £ z a at the acoustic scale, but note that we plot ^ n t/b 2 
and £,sz/b: the linear bias factor of luminous red galaxies is 
b ^ 1.5 - 2.0 lEjlfuifHIfrl. Assuming that galaxies have 
no velocity bias v g = v, the linear bias factor at high red- 
shift is b(z) — 1 = (60 — 1)/D(z), sufficient for £j llt to show 
its structure, when combined with £ zz , £g z , and £ z $, yet the 
plot without b captures the main structure of the correlation 
functions, since the linear bias factor is still of order unity. 

Note that since the source galaxies are at the same redshift 
in Fig. [Hi to Fig.QJi, the cosine angle of the comoving separa- 
tion is 7 = (X2 — Xi)/ r 3 = 0, i- e -> tne redshift-space correla- 
tion function (the sum of the solid, dotted, and dashed lines) is 
6,-dtot = = Co -(1/2) 6 + (3/8) £4, dif- 

ferent from the angle-av erag ed (monopole) correlation £0 of- 
ten used in the literature |j65ll66ll . Figure|2]illustrates the mul- 
tipole components of the redshift-space correlation function at 
z = 0.35. The monopole {solid) is identical to £; nt in shape 
but differs in normalization by £0 = £j nt • cq with the multipole 
coefficient cq in Eq. (T52t . The quadrupole^ {dotted) is nega- 
tive by the sign convention and the hexadecapole £4 {dashed) 
is positive in the figure, while the monopole £0 changes its 
sign as £; nt changes at > 128/i _1 Mpc (see the inset). 



As noted before, the spherical Bessel functions j2{kr^) and 
ji{kr^) in £2 and £4 peak at scales different from the typical 
scale k ~ l/r$ for £o> and thus the acoustic structure seen in 
£0 is smoothed out in £2 and £4. 

In practice, galaxy redshift surveys have a narrow but 
nonzero radial window function and galaxy pairs in the same 
redshift bin often have the line-of-sight separation compara- 
ble to the transverse separation, i.e., 7 7^ 0. The angular de- 
pendence of the redshift-space correlation function, therefore, 
complicates the interpretation of its measurements, which are 
further plagued by low signal-to-noise ratios in estimates of 
£2 and £4. While the monopole £0 can be used to ease the 
theoretical and/or observational challenge, full analysis of the 
anisotropic structure of £ z -dist could in principle bring more 
information than £0 measurements (see 11241 12511 for recent 
analysis). We analyze the full anisotropic structure of the ob- 
served correlation function £ b s below. On small scales, virial 
motions of galaxies result in additional anisotropic structure 
in £ z _dist, known as the Finger-of-God (FoG) effect. Note 
that since this effect involves galaxy motions in nonlinear ob- 
jects, it is not considered in our calculation and linear theory 
provides an inaccurate description of the FoG effect: while a 
simple dispersion model lloTll is often adopted to extract addi- 
tional information contained in the anisotropic structure, it is 
demonstrated ll68ll that this model leads to an unphysical dis- 
tribution of pairwise velocities. However, this difficulty could 
be tackled by recent approach based on modeling nonlinear 
galaxy bias in redshift-space ||69[ F70L F7TT1 . 

The short dashed and short dot-dashed lines in Fig.[T]show 
the correlation functions of the magnification bias £ m b/(5p — 
2) 2 and the evolution bias t; evo /E 2 , respectively. The mag- 
nification bias £ m b is typically smaller than £i nt by the ratio 
of the transverse correlation scale l/k± to the Hubble dis- 
tance 1 /H, and the magnification bias factor is of order unity, 
(5p — 2) = —1.0 ~ 2.0 B28I1 for galaxies and quasars, while 
it can be further suppressed by the source effect canceling the 
volume effect (hp — 2) ~ 0. Since £ m b in Eq. d45l l is pro- 
portional to the projected correlation function w p , its overall 
shape is similar to £i„t but £ m b is positive due to projection 
of £int- As the source population is located at higher red- 
shift, longer line-of-sight distance increases the gravitational 
lensing effect and £ m b increases in redshift, as opposed to 
£i n t oc D 2 oc 1/(1 + z) 2 decreasing in redshift. For exam- 
ple, £ m b at z — 2.5 in Fig.QJ can grow up to a few percent of 
£i nt at the acoustic scale l3ll[32Tl . 

The evolution bias <5 CV o is often ignored in the litera- 
ture compared to the redshift-space distortion bias <5 Z , since 
S cvo oc e ~ V and V <C S z , However, the evolu- 
tion bias can be significantly enhanced when the mean num- 
ber density of sources changes rapidly in redshift. To es- 
timate the evolution boost factor E(z) in Eq. d43b . we as- 
sume the standard functional form of a redshift distribution 
n(z)dz oc z a exp [— (z/zq)^] dz and take two source popu- 
lations as illustrative examples: galaxies and quasars charac- 
terized by (a, P, zq) = (4, 4, 0.4) and (3, 13, 2), respectively 
jl2ll28ll72ll . While the bright samples of luminous red galax- 
ies in the BOSS will have a redshift distribution flatter than the 
assumed here, the faint samples with larger number density 
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FIG. 3: Redshift distribution of source populations and boost factor 
E(z) of the evolution bias. Assuming the standard functional form 
n(z)dz oc z a exp [— (z/«o)H dz, the redshift distribution of galax- 
ies (dashed) and quasars (solid) are shown in the upper panel with 
(a, (3, zo) = (4, 4 ; 0.4) for galaxies and (3, 13, 2) for quasars, re- 
spectively [12, 28, 72]. The bottom panel shows the evolution bias 
boost factor computed by using Eq. 043 1 . The number density of 
sources changes exponentially beyond the mean redshift and the evo- 
lution bias is substantially enhanced in proportion to E(z). 



and volume will have a non-flat redshift distribution lf73lp74ll . 
The clustering of Lyman-a forests at z — 2.5 will be mea- 
sured by the spectrum of quasars at z > 2.5, not by quasars 
themselves at z = 2.5. However, we simply assume that £ b s 
is measured from the galaxy samples at z = 0.35 and z — 0.6, 
and from the quasar samples at z = 2.5. 

The upper panel of Fig. [3] illustrates the redshift distribu- 
tion of the galaxy (dashed) and quasar (solid) samples, with 
its peak at z = 0.4 and 1.8, and the bottom panel shows the 
evolution boost factor E(z) of each sample. For the assumed 
redshift distribution, the evolution boost factor is typically a 
factor ~ 10, and it vanishes at the peak redshift. However, 
a sharp decline in the mean number density of source pop- 
ulations beyond the mean redshift makes the evolution bias 
#evo sensitive to the change in observed redshift z Q b s due the 
generalized Sachs-Wolfe effect, and E(z) can be further en- 
hanced by another factor of ten. With significant boost of 
E 2 (z) ~ 100 - 10000, the correlation £ ovo of the evolution 
bias should be given a proper consideration, especially when 
the mean number of the source population changes rapidly. 
Note that the evolution bias £ cvo /E 2 (short dot-dashed) in 
Fig-Hfc is comparable to the magnification bias £ m b/(5f> — 2) 2 
(short dashed) and is larger at the acoustic peak scale, and the 
evolution boost factor E(z) can be significantly larger than 
the magnification bias factor (5p — 2). Therefore, it is of 
particular importance to select samples of source populations 
that have relatively flat redshift distribution in number density 



(E ~ 10), and to limit the redshift range of measurements be- 
low the peak redshift. The short and the long dot-dashed lines 
in Fig. |2] show the multipole components £ l cvo of the evolution 
bias, defined as £ cvo = J2i=o,2 Ccvo( r 3)fl(7)- B ° th compo- 
nents are smooth and change little over r-j = 1 — 200/i _1 Mpc. 
With£ ovo oc H 2 a 2 D 2 oc 1/(1 + z), it decreases slowly in red- 
shift, and the nonlinear effect is relatively small compared to 
£i n t, since less weight is given to short wavelength modes. 

Now recall that there are four terms of S z in Eq. (|37T i that are 
ignored in our calculation, and they are comparable to e ~ V, 
albeit smaller than the dominant term in 5 Z ~ — tj^t^t- 
While there is no additional boost factor like E(z) in <5 evo , 
their contributions to £ b s are typically of order £ cvo /-E 2 and 
are as large as £ m b at z — 0.35 in Fig.QJz. Though £ m b is 
larger at higher redshift, their impact on £ b s also increases in 
redshift: approximately at the sub-percent level for each con- 
tribution at z = 2.5. This level of accuracy would be appro- 
priate given the statistical errors present in current samples, 
but further calculations of the ignored terms may be needed in 
future surveys. 

In Fig. Q~|i, we consider the correlation functions of two 
source populations, separately located at z\ = 0.35 and zi = 
0.6 as a function of 2D projected separation R = r(x)A6*. 
Note that given a large line-of-sight separation ~ 600/i~ 1 Mpc 
between Z\ and z% , all the correlation functions that depend 
on 3D comoving separation are nearly zero, i.e., £j nt ~ 
£ S z ^ £,zS ^ £zz ^ £cvo ~ 0. The two non-vanishing 
contributions in Fig. Q]i are the auto-correlation of the mag- 
nification bias £ m b/(5pi — 2)(5p2 — 2) (short dashed) and 
the cross-correlation of the intrinsic source fluctuation and the 
magnification bias £5 m b/&i(5p2 — 2) (long dot-dashed) that 
depend on the projected separation, rather than 3D separation 
itself. Note that the cross-correlation in Fig. [TJi to Fig. Q}i 
is identically zero: £g m b = at z\ — 22 with the Limber 
approximation we adopted here, but it is in general smaller 
than £ m b unless z\ 7^ Z2 (see 13 111 for a somewhat different 
derivation). Note that while both £ m b and £s m b in Eqs. ( l45b 
and d55l ) depend on the projected separation via w p , £5 m b has 
additional linear dependence on the comoving line-of-sight 
separation X2 ~ Xi = &X< and it increases with A%, as op- 
posed to £ m b with little dependence on A%. With the large 
A% ~ 600/i _1 Mpc in Fig. Q~|i, £5 m b is substantially larger 
than£ mb . 

Figure|4]examines the anisotropic structure of the observed 
correlation function £ b s , evaluated at z = 0.35. The a>axis 
represents the transverse separation R = r(%)A0 and the y- 
axis represents the line-of-sight separation Ax = X2 — Xi 
with fixed z = (zi + z 2 )/2 = 0.35, x = (xi + Xa)/2 = 
980/i _1 Mpc, andx(^i) < x{ z 2)- The color maps are linearly 
proportional to the value of the correlation function £ plotted 
in each panel below the adopted threshold ( = 4x 10~ 4 , and 
they are logarithmically proportional to £ above the thresh- 
old. The solid contours are also logarithmically spaced with 
increasing thickness at £ > 1.5 x 10~ 3 to emphasize the 
structure shown as the color maps, and their contour values 
are labeled in the color bar. The thickest solid contours sep- 
arate the regions with £ > from those with £ < 0, and the 
dot-dashed and dotted curves represent the contours with £ = 
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FIG. 4: (color online) Anisotropic structure of the observed correlation function £ b s at z = 0.35. We plot the correlation function £i nt of the 
intrinsic source fluctuations in Panel (a), and show the individual effects of the magnification bias fj n t + £mt> + £i mt> (Panel b), the evolution 
bias ^i n t + £ C vo (Panel c), and the redshift-space distortion bias £ z _dist = £int + £zz + £z6 + £sz (Panel d) on the anisotropic structure. The 
observed correlation function £ bs, the sum of all the correlation functions, is shown in Panel (e). We assume that the galaxy bias factor is 
b — 2, the magnification bias factor is (bp — 2) = 2, and the evolution boost factor is E — 100. The color maps are linearly proportional to 
the value of correlation function £ in each panel at £ < 4 x 10~ 4 and are logarithmically proportional to £ at £ > 4 x 10 -4 . The solid contours 
are logarithmically spaced at £ > 1.5 x 10 -3 and their thickness increases with the value of £ (contour values are labeled in the color bar), 
while the thickest solid curves represent the £ = contour. The dot-dashed and dotted curves represent the contours with £ = —4.5 x 10~ 4 
and —9.0 x 10~ 4 , respectively. The acoustic scale in £j nt (ra = 106/i _1 Mpc) is shown as dashed lines in each panel for reference. Note that 
two regions are underrepresented in the color maps, as the region with £ > 1.0 is highly concentrated at r 3 <C 20/i~ 1 Mpc and there is no 
distinctive feature in the anisotropic structure around the region with £ < 0. 



—4.5 x 10~ 4 and —9.0 x 10~ 4 , respectively. For reference, 
we also plot the acoustic scale in £i nt (r$ = 106/i _1 Mpc) as 
dashed lines in each panel. In Fig. Hp,, we plot the correlation 
function £i nt of the intrinsic source fluctuations, assuming the 
galaxy bias factor b = 2. As the rings of the concentric con- 
tours show, £j nt is spherically symmetric and depends only 
on 3D separation r^. The acoustic peak shows its structure 
as a circular ring at r% — 106/i _1 Mpc (dashed) and beyond 
r3 ~ 128/i _1 Mpc £i n t becomes negative without further dis- 
tinctive feature in its structure. 

The gravitational lensing effects, the magnification bias 
£ m b and its cross-correlation £s mb> break the spherical sym- 
metry in £i nt , and its impact on the anisotropic structure 
is shown in Fig. |4p, assuming the magnification bias factor 
(bp — 2) = 2. £ m b depends on the line-of-sight separation A\ 
only through xi an d Xi m Eq. ( 1451 ). and Ax is small com- 
pared to the line-of-sight distance, i.e., Ax <C Xi — X2- Thus 
Cmb is virtually independent of Ax and is just a function of 



transverse separation R, decreasing with increasing R. As is 
seen in Fig.QJi, £ m b is in general orders of magnitude smaller 
than £int at z = 0.35, but £i nt becomes negative at large r$ 
and smaller than £ m b, e.g., £; nt = —8.4 x 10~ 4 < £ m b = 
4.8 x 10~ 6 at R = 10.0h -1 Mpc and Ax = 140/i" 1 Mpc 
(ra = 140.3/i _1 Mpc). Since £ m b changes slowly with R, 
the demarcation curve between the regions with £ m b > £i n t 
and £ m b < £i n t roughly corresponds to the £j nt = contour 
(thickest solid) in Fig. [Ji. However, since |£j nt | > £ m b in 
general except at a narrow strip around the £- mt = contour, 
the impact of £ m b on the anisotropic structure is negligible at 
z = 0.35. Note that the impact of £ m b is substantially en- 
hanced at higher redshift, where longer line-of-sight distance 
results in more fluctuations and the clustering amplitude of 
£ int is lower. 

While both £ m b and £s m b are proportional to w p , (;$ m b de- 
pends on w p itself, rather than the integral of w p along the 
line-of-sight, on which £ m b depends: £g m b becomes nega- 
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tive at large transverse separation R ~ 110/i _1 Mpc as Cint 
becomes negative at large 3D separation r 3 . Note that for 
an observable angular separation A8, the transverse separa- 
tion R = r(x)A6* in the figure is slightly different from 
r(xi)A0 for w p in Eq. $55[ , and this difference tilts otherwise 
a vertical line with £5 m b = at R ~ 110/i _1 Mpc toward 
larger i? at large Ax (6 m b = at R ~ 128/i -1 Mpc and 
A% = 140/i _1 Mpc), because Ax = X2 — Xi and £1 < %2 
with fixed z = (z\ + Z2)/2 = 0.35. Furthermore, since Ca m b 
linearly increases with Ax (hence it vanishes at Ax = 0), 
the absolute value of C5 m b is larger than Cmb except at the 
regions at Ax = and around the nearly vertical strip with 
£,s mb = 0, and it is also comparable to Cint at large Ax and 
small R, e.g., & rab = 4.3 x 1(T 4 at R = 10fc -1 Mpc and 
Ax = 140/i _1 Mpc. Therefore, when Cint is combined with 
Cmb and £g m b as shown in Fig. [4p, Cmb has little impact but 
C<5 mb distorts the symmetric contours of Cint at large Ax and 
small R. At higher redshift, the amplitude of C<5 m b decreases 
with that of £; nt , and Cmb becomes a more dominant contribu- 
tion than C<5 m b- Note that the amplitude at the acoustic scale 
(dashed) is not significantly altered by the gravitational lens- 
ing effects, even along the line-of-sight direction at z = 0.35. 

The impact of the evolution bias £ cvo is shown in Fig. 
where we assume the evolution boost factor E = 100. While 
£g VO and Ce™ change slowly with separation, they have the 
opposite sign as shown in Fig. [2] Therefore, they tend to can- 
cel out along the line-of-sight direction (7 = 1), reducing the 
amplitude of £ cvo , and the largest contribution of £ cvo arises 
along the transverse direction (7 = 0), where the absolute 
values of the monopole Ce VO and the quadrupole Q vo add up. 
For example, £ cvo = 1.6 x 10~ 3 > Cint = -8.4 x 10~ 4 
at R = 140/i -1 Mpc and Ax = 10/i -1 Mpc. As noted in 
Fig. [TJi, the amplitude of £ ovo with E = 100 is smaller than 
£i nt at r3 < 50/i _1 Mpc and its impact is appreciable only at 
t\3 > 100/i _1 Mpc along the transverse direction. As the an- 
gular separation becomes small with fixed 3D separation, the 
impact of £ GVO decreases, because the second order Legendre 
polynomial is a monotonic function of angle. Note that com- 
pared to Cint, £ cvo changes slowly in redshift and its impact 
is larger at higher redshift for a fixed E. The overall shape 
of the acoustic scale (dashed) also remains unaffected by the 
evolution bias. 

Figure |4}i? examines the redshift-space correlation function 
Cz-dist = Cint + Czz + C<5z + Cz<5- Though its angular struc- 
ture is similar to Ccvo, it differs in two aspects: Cz-dist has 
the additional hexadecapole £4, and the quadrupole £2 be- 
comes dominant over the monopole Co at ~ 50/i _1 Mpc 
smaller than r$ ~ 200/i _1 Mpc, where Ce VO ~ Ccvo' shown 
in Fig. |2] Therefore, the angular structure changes more dra- 
matically than that seen in Fig. 0J;, Note that while the fourth 
order Legendre polynomial is not a monotonic function of an- 
gle, the hexadecapole £4 is generally smaller than Co and C2, 
and its contribution is minor. The contours exhibits the well 
known Kaiser effect [20] that coherent infall toward the over- 
dense regions squashes the clustering amplitude and the un- 
derdense regions inflate along the line-of-sight. The large re- 
gion with negative values at Ax — 50— 100/i —1 Mpc and R < 
60/i _1 Mpc is the characteristic feature of this effect, largely 
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FIG. 5: Dissection of the observed galaxy power spectrum P b s - 
Solid and dashed lines represent power spectra of the intrinsic source 
fluctuations Pint/6 2 and the evolution bias P cvo /E 2 . Power spec- 
trum of the magnification bias P m b/(5p — 2) 2 is shown as dotted 
lines for a survey window of width 200ft -1 Mpc (see the text). Note 
that we omit power spectra of the redshift-space distortion bias P zz , 
Ps z , and P z< 5, since they have the same shape as P nt for the line-of- 
sight component of the wavenumber k z — k up to numerical factors 
of order unity. Thin and thick lines represent power spectra com- 
puted by using the linear and nonlinear matter power spectrum. Left 
panel: power spectra are computed at two different redshifts, only 
the power spectrum of the magnification bias increases at higher z, 
while P n t and P cvo decrease. Right panel: cross-power spectra are 
computed for two galaxy populations at Z\ = 0.35 and Z2 = 0.6, 
and dot-dashed lines show the cross-power spectrum of the intrinsic 
source fluctuation at z\ and the magnification bias of source galax- 
ies at Z2- Large line-of-sight separation ~ 600/i~ Mpc corresponds 
to k z — 0.002/iMpc -1 and it has little impact on power spectrum 
(k ~ k± > k z ). 



due to the negative quadrupole C2, and this structure has been 
recently measured with high signal-to-noise ratio 1I24I l25ll . 
Even along the line-of-sight direction, the monopole Co briefly 
takes over the negative quadrupole around the acoustic scale, 
because the clustering amplitude is enhanced. Furthermore, 
while the clustering amplitude increases with angle at the 
acoustic scale (dashed), its structure manifests itself as ridges 

Figure HJ; puts together our discussion, plotting the ob- 
served correlation function Cobs at z — 0.35, the sum of all 
the correlation functions shown in each panel. The redshift- 
space distortion affects the anisotropic structure by far the 
most among the other effects considered here. The gravita- 
tional lensing effects, mostly from Ca mb at low redshift, be- 
come important only at a small transverse but large line-of- 
sight separation. With the same dependence on quadrupole, 
Covo follows the similar angular pattern of Cz-dist> boosting 
their contributions along the transverse direction, but its sole 
impact shows up at R > 110/i _1 Mpc due to the lower am- 
plitude. The clustering amplitude of Cobs at the acoustic scale 
(dashed) is no longer constant, nor a monotonic function of 
angle. The peak location, imprinted in Cmt with local enhance- 
ment in clustering amplitude, remains largely unaffected as 
the gravitational lensing and generalized Sachs-Wolfe effects 
distorts the anisotropic structure. However, note that it is be- 
yond our scope of the current investigation to what accuracy 
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the acoustic pea k remains unaffected by these effects (see, 
e.g., JH HI III M, IH M, HI for recent work on the 
robustness of the baryonic acoustic peak). 

Figure [5] shows the equivalent dissection of the contribu- 
tions of the gravitational lensing and the generalized Sachs- 
Wolfe effects to the observed galaxy power spectrum P Q b s 
in Fourier space. The solid, dashed, and dotted lines repre- 
sent Pint/b 2 , P cvo /E 2 , and P m b/(5p — 2) 2 , respectively. The 
power spectra are also computed by using the linear {thin) 
and the nonlinear {thick) matter power spectrum as in Fig.Q] 
In Fig.|5]2, the source galaxies are assumed to be at the same 
redshift (z = z\ = Z2), and the two sets of lines show the 
power spectra at z — 0.35 and z = 2.5, which decrease in 
redshift except that P m b increases as we have seen £ m b in- 
crease in redshift. The power spectrum of the intrinsic source 
fluctuations P m t/b 2 (solid) exhibits two characteristic scales 
in its structure: a series of the acoustic oscillations starting 
at k ~ 0.085ft.Mpc -1 , and the turnover in the overall shape 
at k ~ Cl.015ft.Mpc~ 1 imprinted by the horizon size at the 
matter-radiation equality z — 3300. 

For simplicity, the wavenumber is set equal to the line-of- 
sight direction k = k z for plotting P cvo and to the transverse 
direction k = k± for plotting P m b. The power spectrum of the 
evolution bias P cvo /E 2 (dashed) is typically many orders- 
of-magnitude smaller than Pmt/b 2 at k > 0.03ftMpc _1 in 
Fourier space, but its contribution can be at the few percent 
level of Pi n t at the acoustic scale and comparable to P; nt at the 
matter-radiation equality scale, with the evolution boost fac- 
tor E 2 ~ 10000. Note that the power spectra of the redshift- 
distortion bias P zz , P$ z , and P z g are omitted in the figure, be- 
cause they have the same shape as p nt up to numerical factors 
of order unity when k = k z . 

To plot the power spectrum of the magnification bias 
P m b/(5p— 2) 2 (dotted), we replace (27r)(5 D (fc z ) by aflat win- 
dow function of width 200/i _1 Mpc, typical value in redshift 
surveys, hence the dotted line de facto delineates the angular 
power spectrum of the magnification bias (5p— 2) 2 C*J? kj _ r ,-s 

with a dimensional coefficient r 2 (x) X (200/i -1 Mpc) (see 
Appendix [A}. While the magnification bias is negligible at 
z = 0.35, its effect increases with larger line-of-sight distance 
at higher redshift, amounting to a few percent at the acous- 
tic scale and larger at the matter-radiation equality scale at 
z = 2.5. However, note that even with relatively large con- 
tributions to P bs, the shift in the peak positions can be at the 
sub-percent level or smaller ll3lll32tl . 

Figure|5l) plots the cross power spectra of two source popu- 
lations at z\ = 0.35 and Z2 — 0.6. As opposed to the correla- 
tion functions shown in Fig.[T]i, all the power spectra that de- 
pend on 3D wavenumber remains virtually unaffected by the 
large line-of-sight separation ~ 600/i _1 Mpc, because it cor- 
responds to very small wavenumber k z ~ 0.002/iMpc~ 1 and 
k ~ k± ^s> k z . The dot-dashed lines show the cross power 
spectrum of the intrinsic source fluctuation and the magnifi- 
cation bias Ps m h- Since it is proportional to the line-of-sight 
separation, its contribution is larger than P m b in Fig.|5f), but it 
is absent in Fig.[5]2. 

The anisotropic structure of the observed power spectrum 
P bs has been well studied with main focus on the effect of 



the redshift-space distortion bias I20I1 . and the redshift-space 
power spectrum P z -dist in Eq. d53l ) has the multipole compo- 
nents that are identical in shape but only differ in normaliza- 
tion. The evolution bias results in the similar angular pattern: 
two multipole components that share its shape with p nt with 
different normalization. However, note that since the magni- 
fication bias and its cross term are intrinsically 2D quantities, 
their impact on the anisotropic structure of P Q b s is small, even 
with realistic survey window functions P 55ll . 



V. DISCUSSION 

Galaxy two-point statistics, correlation function in real 
space and power spectrum in Fourier space, have been ex- 
tensively used in cosmology to characterize the underlying 
matter fluctuations. We have presented a coherent theoreti- 
cal framework based on the linearized Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric for computing the gravita- 
tional lensing and the generalized Sachs-Wolfe effects. Within 
this framework, the metric perturbations are sourced by the 
underlying matter fluctuations, and they naturally give rise to 
perturbations in the observable redshift of source galaxies and 
their angular position on the sky. The time component of the 
photon geodesic equations can be used to show the former, 
the generalized Sachs-Wolfe effect [40] that generalizes the 
standard redshift-space distortion by peculiar velocities in a 
cosmological context, including the Sachs-Wolfe and the in- 
tegrated Sachs-Wolfe effects. The spatial components of the 
photon geodesic equations can be used to derive the latter, 
the gravitational lensing effect that includes the weak lensing 
distortion, magnification, and time delay effects. This uni- 
fied treatment provides a complete description of the relation 
between these seemingly different effects and the underlying 
matter fluctuations. 

Furthermore, it becomes transparent in this treatment how 
the gravitational lensing and the generalized Sachs-Wolfe ef- 
fects affect the observed fluctuation field of source galaxies. 
To the linear order in perturbations, we have computed all 
the additional contributions to the intrinsic source fluctua- 
tion, arising from the gravitational lensing and the generalized 
Sachs-Wolfe effects. We can gain more insight on the impact 
of these effects by separating them as two physically distinct 
origins: the volume and the source effects. The former effect 
that involves the change of volume is independent of source 
galaxy populations and hence regardless thereof the volume 
effect is always present in galaxy two-point statistics. By con- 
traries, the latter effect depends on the intrinsic properties of 
source galaxy populations and may vanish for a certain popu- 
lation. All of these contributions to the intrinsic source fluc- 
tuations result in numerous additional auto and cross terms in 
the observed galaxy two-point statistics, and therefore proper 
account should be taken into these additional terms in inter- 
preting measurements of galaxy two-point statistics from up- 
coming dark energy surveys. 

With the complete list of the contributions of the gravita- 
tional lensing and the generalized Sachs-Wolfe effects, sepa- 
rated as two physically distinct origins, we have identified sev- 
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eral contributions in the volume effect and one contribution in 
the source effect, which are ignored in the standard treatment: 
the evolution bias in the source effect arises from the gener- 
alized Sachs-Wolfe effect, when the mean number density of 
sources changes rapidly in redshift, and its impact on the ob- 
served galaxy two-point statistics can be substantially larger 
than that of the gravitational lensing magnification bias. The 
ignored contributions in the volume effect are typically of or- 
der peculiar velocities and hence they are subdominant, com- 
pared to the standard redshift-space distortion effect. How- 
ever, their impact is comparable to the magnification bias at 
low redshift. While the cross term of the magnification bias 
and the intrinsic source fluctuation is more important at low 
redshift than the contribution of the magnification bias itself in 
the gravitational lensing effect, further calculations of the ad- 
ditional contributions associated with the volume effect may 
be needed, if higher accuracy of theoretical modeling is re- 
quired from observation. 

We have investigated the impact of the additional contribu- 
tions to the anisotropic structure of the observed galaxy two- 
point statistics, after simplifying some of the contributions 
to the intrinsic source fluctuations. The redshift-space dis- 
tortion affects the observed galaxy two-point statistics most, 
imprinting its well-known feature in the anisotropic struc- 
ture ll20l l46l l47ll . The gravitational lensing effect is small 
but non-negligible at a percent level, particularly along the 
line-of-sight separation and at high redshift, since their con- 
tribution increases with longer line-of-sight distance to the 
source galaxies and the clustering amplitude of the intrinsic 
source fluctuations decreases in redshift. The evolution bias 
has an angular pattern similar to the redshift-space distor- 
tion, but its impact becomes appreciable, only at fairly large 
transverse separation. While it is challenging to analyze the 
observed anisotropic structure of galaxy two-point statistics, 
its full analysis from upcoming dark energy surveys can pro- 
vide a great opportunity to separately identify each contribu- 
tion from the gravitational lensing and the generalized Sachs- 
Wolfe effects, increasing the leverage to understand the un- 
derlying physical mechanism. 

However, we note that constraining the underlying cosmo- 
logical model will require not only accurate theoretical predic- 
tions, but also model fitting to measurements, which results 
in further distortion in galaxy two-point statistics, known as 
Alcock-Paczyhski effect [81]. Furthermore, our current in- 
vestigation has focused on the linear theory predictions and 
its additional contributions: nonlinearity and scale-dependent 
galaxy bias can affect our results, though its impact is ex- 
pected to be less than at the percent level around the acoustic 
scale (see, e.g., JB M IH). However, additional leverage 
can be gained by modeling scale-dependent galaxy bias on 
nonlinear scales [821. 
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APPENDIX A: 2D AND 3D STATISTICS 

Here we derive the relation between 2D and 3D fluctua- 
tions and their two-point statistics. Consider a fluctuation field 
<5 2D (n; z s ) on the sky from a source population at z s . In gen- 
eral, it can be expressed in terms of the convolution of a win- 
dow function ^(x) and its 3D fluctuation <5 3D (x) 



S 2U (h;z s 



dxW( Xs -x)6 6U Hx)n,X;z}. (Al) 



When the window function is appreciable only around z s 
representing a narrow selection function in redshift surveys, 
<5 2D (n; z s ) ~ <5 3D (x; z s ) with its functional dependence x = 
[r(x s )n, Xs]- However, contributions to <5 2D (n; z s ) can come 
from the fluctuations <5 3D (x; z) at z < z s and <5 2D (n; z s ) may 
be substantially different from <5 3D (x; z s ), when the window 
function is broad. For example, the convergence field fc(n; z s ) 
in Eq. d2Tb has the window function 



w K ( Xs -x) 



r(Xs ~ X) r( X ) 
a(x) r (Xa) 



(A2) 



which peaks roughly at a half of r(x s )- 

In a sufficiently small patch of the sky, the Fourier mode of 

n; z s ) is 



d 2 he- llil 5 2D (h;z s ) 

r 2 (X) 

dkz e ik,X £3D 
27T 



(A3) 



k z , — — ; 



and its (angular) power spectrum is 



Ci(zi 1 z 2 ) 



d 2 r 

dXa I '/\/, 
dk 



(Sf-( Zl )Sr(z 2 )) 

W(xi-Xa)W( X 2-Xb) 



(A4) 



r (Xa) 2 

k z i k± = 



I 



d W{xi-x)W(x2-x) p 
r 2 (x) 



r(Xa 
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I 



r(Xa) 
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The last equality is obtained by adopting the Limber 
approximation, in which fluctuations along the line-of-sight 
are smoothed out and only long wavelength modes (k z ~ 0) 
can contribute to the integral l54l l83ll . With the Limber 
approximation, the angular correlation function is 

w(A6;zi,zz) = (,5 2D (ni;zi) <5 2D (n 2 ;z 2 )) (A5) 

= / d X W(xi-x)W( X 2-x) 
Jo 

x w p [r(x)A0; z] , 



In practice, S D (k £ ) need to be replaced by a survey 
window function 115 511 . but note that it is crucial to assume 
the independence of source redshift, when computing the 
power spectrum. We also note that the Limber approximation 
breaks down when the radial window function of a survey 
is narrow compared to the correlation length scale (see, e.g., 
tfS4l HH]). However, the use of the Limber approximation 
is readily justified in galaxy surveys, in which the radial 
window function has width of Az ~ 0.1 — 0.2, corresponding 
to several hundred Mpc. 



where the projected correlation function is 



APPENDIX B: GAUGE-INVARIANT FORM OF OBSERVED 
NUMBER DENSITY 



w p [R; z] = 



k dk 



= , R 2 +r 2 :z 



(A6) 



2tt 



P(k;z)J (kR). 



The angular correlation function and power spectrum in 
Eqs. ( |45T > and d47l i can be readily obtained by substituting the 
window function W K ( X ) for the convergence in Eq. (IA21 > with 
W( X ) in Eqs. iA4i and jA5\ . The cross correlation function 
and power spectrum in Eqs. (IB3T l and (TS6b can be computed 
in a similar manner, since W(x s — x) — ^ D (Xs ~ x) gives 
<5 2D (n;z s ) = 5 3D (x;z s ). 

In Sec. |IV] we associated the angular power spectrum C" K 
to a 3D power spectrum to compare its impact with other 3D 
power spectra. A 3D fluctuation field can be constructed from 
<5 2D (fi;* s )by 



S(k;z s ) = J d 3 xe-^*S 2D (h;z s ) (A7) 

= J d Xs r 2 ( X s)e- ik *x»5™ k±rs (z s ) 
= (2n)5 D (k z )r 2 (xs)S™ k±rs (z s ). 



We assumed Sf° is independent of z s in the last equal- 
ity. For high redshift source populations, this approximation 
is accurate, since the growth of the comoving angular diame- 
ter distance flattens at high z and it becomes nearly constant. 
Within this approximation r(xi) = KX2) — r (x)> me 3D 
power spectrum is anisotropic and it is related to the angular 
power spectrum by 



P(k z ,kj_;zi,z 2 ) 



d 3 k' 



(«(k;zi)a*(k';* 2 )} (A8) 



2tt) 3 

(2tt) 6 D (k z )r 2 (x)C l=k±rm (z) 



Here we provide a rigorous derivation of the observed num- 
ber density in Sec. lIII Al For simplicity, we assume a flat uni- 
verse. 

In the observer's frame, local coordinates p a are used to 
describe the observed positions of galaxies and their true po- 
sitions are related to the observed positions by the photon 
geodesic x a (X). The total number of observed galaxies can 
be computed by considering a covariant volume integral 18611 : 



/ 



-9 "phy u dS d 



(Bl) 



where the (oriented) hyper-surface element is 



dx a dx b dx c ,1,9,3 
dS d = Eabcd — — — dp dp dp 



dp 1 dp 2 dp 3 



(B2) 



g is the determinant of the space-time metric, n p h y is the phys- 
ical number density of galaxies, and t a bcd = £[ a bcd\ is th e 
Levi-Civita tensor density. We take the Newtonian gauge vari- 
ables (z, 9, 4>) as the observed local coordinates p a . For nota- 
tional simplicity, tilde is used to represent observed quantities 
(e.g., z = z bs)- By imposing the number conservation, the 
observed number density is then related to the total number of 
galaxies by 



N, 



gal 



Y 2 (S) 

n ^4 saxOdzdOdd) 
H{z) 



(B3) 



In a homogeneous universe, the local coordinates are identi- 
cal to the true coordinates (z, 9, (j>) = (z, 9, <fi), and the photon 
geodesic is simply ir a (A) = (y, y e Q ), where we choose the 
normalization of the affine parameter as A = (a/v)y. Noting 
that the four velocity of a comoving observer is u a = (1/a, 0), 
equation (IB 11 1 can be readily solved as 
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Nr. 



4 n phy dx a dx@ dx 1 - ~ 



a n phy 



a dz 

x 2 W 



H(z) 



86 

sin 9 dzd9d(j). 



(B4) 



Therefore, we recover the standard relation for iVg a i 
and h — a 3 nphy = n{z, 9, <p). 

In an inhomogeneous universe, the photon geodesic devi- 
ates from the null path and the local coordinates are differ- 
ent from the true coordinates. Perturbations to the photon 
geodesic in an inhomogeneous universe can be computed by 
integrating the null vector k a (X) = dx a /dX, 



1 fdx a dxPdx~<\ (dz 89 



, (B8) 



where the subscripts denote the order in perturbations, to 
which quantities in the bracket need to be computed. When 
combined together, equation {B6} is 



u dbd = — — ; sin 6 



a H{z) 
1 + Se a e a 



dz 89 
8S + 89 



— + — = + — £ -ip + v a e 



(B9) 



x a (\) = (y,ye a ) + / dy' (6u, Se a ) . (B5) 
J o 



Note that to the first order in perturbations the integration is 
performed along the null path, ranging from the observer at 
y = to the source galaxies at y. 



With u a 
tion (IbTT i is 



((1 — ip) / a,v a / a), the integrand of equa- 



u d dS d 



1 — ip 8x a 8x@ dx 1 



(I 



82 89 



v a 8x a dx 8x c 



Finally, the determinant in equation (IBU gives y/—g = 
a 4 (1 + ip + 3 <p) and the total number of observed galaxies is 



•Y„,i = / o 3 nphy ^4 siae<ll,W<h> (BIO) 



H(z) 

, dz 89 8(h 
dz d9 d(j)J 1 



(B6) From equation ( |B3t , we obtain our final result, 



The last two terms, proportional to ip and v a , contribute to 
the first order in perturbations and the partial derivatives need 
to be computed, only to the zeroth order. The first term has 
two sources of perturbations from the partial derivatives: 
perturbations in the photon geodesic and the relation between 
the local and true coordinates. The former is non-zero, only 
when the derivative is taken with respect to z, i.e., 



1 fdxPdx~t\ Se a 

Cq/370 



a 



89 



o H{z) ' 



(B7) 



and the latter is 



X 2 {z) H{z) 
X 2 (z) H{z) 



1 + 20- (1 



x de 
z)— - 2 k 
dz 



(Bll) 



which includes the distortion of volume element SV and grav- 
itational lensing magnification. This expression is manifestly 
gauge-invariant and valid on all scales. In special relativity, 
the volume element of the local Lorentz frame of matter is dis- 
torted by 7 = vl — v 2 , and hence is identical to an observer 
at rest, to the first order in perturbations. However, in our case 
the volume element can be measured, only by observing light 
rays of photons, of which time component is related to the 
spatial component, giving rise to the first order distortion in 
volume element. 
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